On the mechanism of trapped surface formation 

D. Valls-Gabaud^'^ and T. Zannias"^ 
Department of Physics, Queen's University, Kingston, Ontario K7L 3N6, Canada. 



^ ABSTRACT 
0^ 

bJ[) By a combination of analytical and numerical methods, the density profile of a 

^ momentarily at rest spherical star is varied, and the corresponding response in the area 

of the spherical shells is monitored. It is shown that the inner apparent horizon (if 
0> it exists) must lie within or at most on the star's surface, while no such restriction is 

found for the outer apparent horizon. Bowever, an apparent horizon lying in the vacuum 
region will always have non vanishing area, as long as the ADM mass of the system is 
non zero. Furthermore for density profiles not decreasing inwards, it appears that all 
lO spherical trapped surfaces lie on a thick spherical shell. Finally for a uniform density 

star a simple criterion is found, relating density and proper radius that guarantees the 
presence or absence of trapped regions. 
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There has been lately considerable renewed effort aiming to discover conditions, upon a given 
initial data set, that permit us to infer the presence or absence of trapped surfaces [1,3,4-6]. Al- 
though most of the cited work is focused on characterizing the state of initial data triggering the 
formation of trapped surfaces, it would also be beneficial to understand the mechanism that leads 
to their formation as well as any other relevant properties. For instance how do spherical trapped 
surfaces appear and how are they distributed around the centre of the star? Where (if it exists) is 
the location of the outer apparent horizon? What is the response of the apparent horizon under 
smooth variations of the parameters describing the initial data? 

It is the pupose of the present note to provide, at least for a special class of configurations, 
answers to the above-mentioned questions. For that we shall consider a sequence of instantaneous 
states describing spherical stars of proper radius R and proper density p just at the onset of their 
gravitational collapse. By examining the proper area of spheres around the centre of the star as a 
function of distance away from the the centre, the trapped region can be identified: it lies between 
succesive maxima and minima of the area function. We vary the density profile (in essence we are 
moving from one member of the sequence to another) and monitor the behaviour of the area func- 
tion. Such "variations" of the density cause the appearance or dissappearance of trapped surfaces 
and forces the inner and outer apparent horizons to "move" around in the star's interior or exterior 
regions. 

Recall that initial data on a time-symmetric slice S satify the hamiltonian constraint, i.e., 

= Wtt^p (1) 

where p stands for the non-negative density. The assumption of spherical symmetry allows the 
introduction of geodesic-type coordinates and thus the line element of S takes the form 

d.s'^ = dr^ B{r) (de^ sin^ Odcf)'^^ < r < oo (2) 

Evaluating the scalar curvature ^TZ of the metric (2), equation (1) reduces to 
2 d'^B 1 fdB\^ 2 _ G 
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We look for solutions admitting regular geometry at r = 0. This demands 

dB 



B(r) 



r=0 — — — 

dr 



= (4) 

r=0 



For later use note that the amount of mass-energy within each sphere of radius r is given by [2] 

(4a) 



1 ^ (^^^ ^ 
~ Jb \ d^ 



and it obeys, by virtue of Eq. (3) to 

^ = 2.pB'l'^4- (4b) 
dr dr 
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Our intention is to investigate the behaviour of the solutions of Eq. (3) subject to (4). A number 
of conclusions can be drawn by inspection of Eq. (3). At first for p = a solution obeying (4) is 
given by B(r) = r^, which corresponds to a fiat-three space. Further if a solution B = B{r) admits 
infiexion points of non zero area (and thus the area function passes through a saddle point) they 
must be either interior points or lie at the star's surface. At such points, one deduces from Eq. (3) 
that the density p and the area B{r) obey the relation Bp = c'^ /SttG. 

Let us now consider a solution extending from the interior to the exterior vacuum region. Absence 
of surface layers at the star's surface requires B(r) to be at r = R, i.e., B(r) and dB/dr should 
be continuous across the surface (note however that d^B/dr^ may exhibit discontinuous behaviour 
at r = R). Since on the other hand in the vacuum Eq. (3) implies d^B/dr^ > 0, therefore if 
the surface of the star is not trapped (i.e., dB/dr > there), then the exterior solution will be a 
monotonically increasing function, that is, trapped surfaces can never develop in the vacuum region. 
In the opposite case, i.e., if the surface is trapped, Eq. (3) implies that only a local minimum of 
B(r) can develop in the vacuum region. This outermost local minimum is the outer apparent 
horizon. Note however that in general it may also lie in the interior. In the case where it lies in 
the vacuum one may naturally ask: can it have zero area? We shall show that as long as the data 
possess a non zero ADM mass, this is not possible. In fact evaluating the right hand side of Eq. 
(4a) on the outer horizon (of zero area), and using Eq. (3) we get : 



which, in view of Eq. (4b) is impossible. This argument also rules out vacuum infiexion points of 
zero area. 

Let us now shift attention to the behaviour of the local maxima of B(r). The first local maximum 
of B(r) (if it exists) marks the location of the inner apparent horizon. Generically it lies interior to 
the surface, but one may ask whether the inner horizon lies at the surface of the star. The answer 
depends on the behaviour of the density profile at the surface. If p(r) is continuous at r = R then 
Eq. (3) implies that the inner horizon cannot be located at the surface. If it exists it must be 
an interior point. In the case where p(r) is discontinuous at r = R, Eq. (3) implies that the left 
derivative {d^B/dr^)j^ and surface density satisfy 




Therefore arbitrary density profiles allow the possibility that the left second derivative be negative 
while simultaneously the right hand side of (5) is positive, implying that the surface of the star 
may become the location of the inner horizon. Notice in that case the matching of B(r) across 
the surface dictates that the second right derivative is given by 



implying that the surface must be an outer apparent horizon as sensed by the exterior geometry 
(this situation will be verified below by an explicit exact equation). In summary, no inner horizon 






(5) 




(6) 
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Figure 1: Dimensionless area /3 as a function of the dimensionless proper distance y = r/R for a 
uniform density distribution. The values of the parameter Tq are indicated. 
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Figure 2: As Figure 1, for gaussian density distributions characterized by Tq and a. The dotted 
vertical lines mark the positions of local extrema. Note the presence of a maximum and a minimum 
(indicative of a trapped surface) in the area above a given threshold in the parameter Tq, for a 
given compactness parameter a. 
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can lie exterior to the star in contradistinction to the outer apparent horizon. 



How many distinct local maxima and minima are admitted by an arbitrary solution of Eqs. (3) 
and (4) matched in a fashion across the surface? The answer to this question is not obvious. 
The non-linear structure of Eq. (3) does not allow us to get further insight into the nature of 
the solutions and thus we have to search for exact solutions. For arbitrary density profiles p(r) 
exact solutions appear difficult to find, therefore we shall use numerical techniques to construct 
solutions of the system. Note that one case where the system is soluble corresponds to uniform 
density stars, i.e., p = const. [7]. Before we comment on this solution, let us rewrite the system in 
a slightly different form, more suitable for numerical computations. It is convenient to introduce 
the dimensionless variables 



y = r/R [i = BlR^ 



where R is the physical (proper) radius of the star, so that Eqs. (3) and (4) take the dimensionless 
form 

2>fB 1 fdB\^ 2 _ 



/3 



y=0 



d£ 
dy 



(9) 



y=0 



where 



T{y) = ^R'p(r 



10) 



In this notation, we shall examine two density profiles corresponding respectively to a uniform 
configuration 



ruiy) 



To 0<y<l 
y>l 



111 



and a gaussian one 



e«(i-!/') <y <1 
y > 1 



'12) 



where Tq and the compactness parameter a characterize the profiles. The solution of Eqs. (8) 
(9) with source given by (11) has the form in the interior (y < 1) 



SttT, 



sm 



SttT, 



1/2 



and 



13) 



In the exterior region Eq. (8) is numerically integrated with boundary conditions dictated by 
the matching of the solution on the surface. In Figure 1 solution curves for various values of 
the dimensionless parameter Tq are shown. It is worth to mention briefiy a few features of the 
solution relevant to our discussion. One may easily verify that as long as Tq < 37r/32, (3(y) in- 
creases monotonically. At precisely this critical value Tq = 37r/32, (3(y) develops the first local 
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Figure 3: Examples of configurations for which several maxima and minima are present in the 
dimensionless area (3. (a) Uniform density configuration, (b) Gaussian density profile. 

maximum taking place at the surface of the star. Note that because of the discontinuity of the 
density across the surface, this critical point is a local minimum according to the exterior geometry. 
Thus for this particular value of Tq there exists only one trapped surface which is simultaneously 
the location of the inner and outer apparent horizons. As Tq keeps increasing, the degeneracy is 
lifted and the inner and outer apparent horizons start "moving" on opposite directions (see also 
Fig. 1). Upon further increase, the inner horizon moves moderately towards the centre while the 
outer horizon moves towards the surface and simultaneously keeps shrinking (i.e. its proper area 
decreases). Finally at Tq = Svr/S the outer horizon has zero proper area and resides exactly on the 
surface of the star. Any subsequent increase of Tq causes the outer horizon to move inwards while 
maintaining its zero area. For even larger values of Tq more critical points start to appear (see 
also Fig. 3a), but such data may be regarded as unphysical since the gravitational field generated 
by them is so strong that it disconnects the star from the asympotically fiat region (for further 
relevant comments on this issue see ref. [8]). In summary, whenever data obey 37r/32 < Tq < 37r/8, 
they admit a distinct inner and outer apparent horizons. 

Solutions of Eqs. (8) and (9) with a gaussian density profile (12) exhibit the same qualitative 
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features as those of the uniform case. In Figure 2 (a-f) we present a few solutions for various 
values of the parameters a and Tq. Note however that in this case there exists data that allow 
both inner and outer apparent horizons of non zero area lying within the interior of the star (Fig. 
2e), unlike the solutions for uniform density stars. From the analysis of the solutions it appears a 
common property worth emphasizing. The system of Eqs. (3) and (4) admits no solutions where 
two or more consecutive minima (excluding the origin) have values different from zero. In turn, 
this property implies that a system must first disconnect itself from the asymptotically fiat region 
before a second minimum appears (see Fig. 3 a-b). We like to think that this property of Eqs. 
(3) and (4) is a generic feature of all solutions, but we have been unable to show this analytically. 
We discuss elsewhere [9] the impact of these constraints on models of compact astrophysical objects. 

Finally it is of interest to know how much of the above described picture is maintained if one deals 
with configurations lacking time symmetry or spherical symmetry. Although one may anticipate a 
similar overall behaviour (at least for spherical systems), we ought to bear in mind that there are 
additional factors entering the problem which may alter some of the features discussed here. We 
hope to come back to this point elsewhere. 
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